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The thermal pairwise entanglement (TE) of the 5=1/2 XY chain in a transverse magnetic field 
is exactly resolved by means of the Jordan- Wigner transformation in the thermodynamic limit 
N^oo. It is found that the TE vanishes at a fixed point with temperature T C ~0.4843J, which is 
independent of the magnetic field. A thermal quantity is proposed to witness the entangled state. 
Furthermore, the TE of the 5=1/2 antiferromagnetic-ferromagnetic (AF-F) Heisenberg chain is 
studied by the transfer-matrix renormalization group method. The TEs of the spins coupled by AF 
and F interactions are found to behave distinctively. The vanishing temperature of the field-induced 
TE of the spins coupled by F interactions is observed dependent on the magnetic field. The results 
are further confirmed and analyzed within a mean-field framework. 
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Quantum entanglement describes intrinsic correlations 
incurred in quantum mechanics. It plays the essential 
role in quantum information [l[, quantum teleportation 
Q , and quantum cryptography [3| . In condensed matter 
physics, it provides a new perspective to understand the 
collective phenomena in many-body systems Q. 



The spin entanglement in quantum spin chains is of 
particular interest. Many types of entanglement at both 
zero and finite temperatures have been extensively stud- 
ied in various spin systems (see reviews in Ref. 0). As the 
finite-temperature entanglement (thermal entanglement, 
TE) can be witnessed theoretically [|| and detected ex- 
perimentally [6( by macroscopic variables, most effects 
have been made to quantify the TE. The critical tem- 
perature (CT) below which the TE survives can even be 
estimated in experiment Q. Interestingly, the numerical 
calculations indicate that the CT of the nearest-neighbor 
TE is a fixed point which is independent of the magnetic 
field in the 5= 1/2 Heisenberg chain Q and two-qubit 
XY spins However, the features at the fixed point are 
unclear and no explanation exists. It is also a question 
whether the magnetic field independence of the CT is a 
universal phenomenon in quantum spin chains or there 
are exceptions. For these questions, in the paper, the 
field dependence of the CT of TE in the 5=1/2 XY and 
AF-F Heisenberg chains are exactly resolved and studied 
by means of the transfer-matrix renormalization group 
(TMRG) method in the thermodynamic limit N^oo, 
respectively. An analysis will also be made within the 
mean-field framework. 



The pairwise entanglement of two 5=1/2 spins at sites 
i and j in the ground state and at finite temperature 
can be achieved from the corresponding reduced density 
matrix pij, which, in the standard basis {|TT), |TI)? lit 
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where Pl=±(l+a z ), P^=i(l-a z ), and a ± =±(a x ±a y ). 
The brackets denote the ground-state and thermody- 
namic average values at zero and finite temperatures, re- 
spectively, and a are Pauli matrices. As the phenomenon 
that is of interest mainly exists in the nearest-neighbor 
TE, we shall concentrate only on Pi^+\ in the following. 

The spin operators can be transformed into spinless 
fcrmions by the Jordan- Wigner (JW) transformation 
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where cj and Cj are the creation and annihilation opera- 



tors of the spinless fermion, respectively. 
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where X i f =(n i n i+ i) (^=4^), 

Yf={rii+i(l - rn)), Zi={c\c i+ i), and XT=l - (m) - 
(rii+i) + (rijni+i). As defined, the concurrence of TE of 
two nearest neighbors is given through 



Cj = Hi - p,2 - fJ>3 - M4j 

Cj = max{0, Ci}, 



(4) 
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where /ij are the square roots of the eigenvalues of 
pi,i+xpi,i+\i where p,\ is the largest, pi.i+i is a trans- 
formed matrix of pi,i+i, i.e., p=a y ®a y p* a y ®a y . Thus, 
Eq. (J4j) is transformed into 
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The concurrence can be calculated from the local density, 
hopping term, and site-site correlations of the fermions. 

As the observed field independence of the CT of TE are 
obtained by either numerical calculations @ or only for 
two qubits [9( , a deep understanding is indeed necessary. 
Therefore, we shall exactly resolve the concurrence of the 
5=1/2 XY chain in a transverse magnetic field within 
the thermodynamic limit to investigate this phenomenon 
analytically. The Hamiltonian of the 5=1/2 XY chain is 
given as 

N N 

H = Y J ^J{sis; +l + h. c .)-hY J S z l , (7) 

where J(>0) is the coupling, and h is the magnetic field. 
As the XY chain with F couplings can be obtained by a 
unitary transformation, which rotates the odd-site spins 
by 7r angle around the z-axis, both the AF and F cases 
give the same results. Here we take J>0 for simplicity. 

By applying the JW and Fourier transformations, the 
Hamiltonian (0 can be diagonalized as 

H = Y^{J ™sk - h)c{c k =Y^Hk) - h]c\c k , (8) 

k k 

and the elements in the reduced density matrix [Eq. l[5])] 
I 



Q = 



x 



where f3'=/3J and h'=h/J. The result of Eq. {THJ) is 
shown in Fig. [TJ where the TEs in different fields vanish 
at a common CT (T c ), which is a fixed point. This com- 
mon CT indicates that the magnetic field cannot retrieve 
the intrinsic TE once it is destroyed by thermal fluctu- 
ations, even though the field changes the TE below the 
CT. 

I 



can be explicitly expressed as 
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where f(k)=l/(e^ E ^~ h ^ + 1) {(3 is the inverse temper- 
ature and the Boltzmann constant is taken as fcs = l) is 
the Fermi distribution function. By the solution of the re- 
tarded Green's function G r (t)=<C c kl {t)ck 2 (t), c£ c£ ^> 
(fciT^) [l(|, the expectation value (c r ki c i k2 Ck 1 Ck 2 ) in Eq. 
(flC)]) is obtained as 

(4A 2 c ^ 2 ) = -f(h)f(k 2 ) (h^k 2 ), (ii) 

and Eq. (I10| is simplified as 

X+ = («,) 2 - Zl (12) 

By substituting Eqs. © and $T2§ into Eq. ©, C { can 
be obtained as 



I 

At T=T C , Ci=0 and Z?=X+Xr. Thus, we can derive 
the equation 

(m) - (n,) 2 = -V2Z t - Zf (14) 

at T c for any fields. We define $(/?, ft,) = (n l )+ v / 2Z i +Z J 2 - 
(rti) 2 , which can be expressed as 
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Thus, <I><0 describes the entangled state, and $ must 
be zero at T c , i.e., $(/3 c ,/i)=0, from which T c can be 
determined. In the absence of magnetic field, (rij)=l/2 



at any temperatures, thus, Zj a h=0 ={l — \/2)/2, and T c 
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FIG. 1: (Color online) Temperature and magnetic field depen- 
dence of the thermal entanglement in the S=l/2 XY chain. 
The entanglements vanish at a common temperature T c . 



can be obtained by solving the equation 
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which indicates that T c is proportional to J, i.e., T c =aJ 
with a~0.4843. At T c , <&(/3 c ,/i) is independent of h, i.e., 
<9$(/? c , h)/dh—0, yielding the equation 
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which is satisfied at T c for any fields. This equation, as 
well as Eq. (fl"4"]) . determine the fixed point completely. 
In Ref. H, the CT for the 5=1/2 two cyclic XY qubits 
is 0.5673J, which is larger than the CT 0.4843J in the 
thermodynamic limit. This is consistent with the result 
in Ref. 8, where the CT of the spin- 1/2 Heisenberg chains 
with iV=5^10 are smaller than that with N=2. 

In terms of the spin operators, the finite pairwise TE 
survives when 
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which indicates that the concurrence of TE is determined 
by the competition between the spin fluctuations and lo- 
cal magnetic moment at finite temperature. As the quan- 
tities in Eq. ^ can be expressed as thermodynamic ob- 
servables as 
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where U=(H) is the internal energy, and M—^^Sf) is 
the total magnetization, the TE can be witnessed by the 
negative thermal quantity 
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FIG. 2: (Color online) Temperature dependence of the ther- 
mal entanglement of (a) C a , and (b) C/ for the 5=1/2 AF-F 
chain at various fields obtained by means of the TMRG. 



which can be measured in experiment. Without the field, 
the magnetization M vanishes, and the TE survives when 



NJ 2 1 



(21) 



which includes a wider parameter region than the suffi- 
cient condition \U\/N J>\/A for the entangled state that 
is proposed for the spin chains with Heisenberg or XY in- 
teractions 0]. In a magnetic field, Eq. ([20]) also implies 
that the witness \U + Mh\/NJ>l/4 proposed in Ref. 5 
can be improved to cover wider parameter region for the 
entangled state. The exact solution not only confirms 
that the CT of the intrinsic TE which survives in the 
absence of field is a fixed point, but also reveals some 
features at the fixed point from the perspectives of the 
local spin competition and macroscopic thermodynamic 
behavior. 

To investigate the possible exceptions of the field inde- 
pendence of the CT, we study the TE in an 5=1/2 AF-F 
alternating Heisenberg chain by means of the TMRG. 
The Hamiltonian of this alternating chain is given by 



h — y^(j a s2j-i 



Jf&2j ■ &2j+l) — h > 5,f 
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(22) 



where J a >0, Jf<0 denote the AF and F couplings, re- 
spectively. J a is taken as the energy scale and Jf/J a =— 1. 
This AF-F chain has a Haldane gap A~0.6J a in the 
ground state [ll|, and the saturation field h s ~l.U a . In 
experiment, this model has been realized and studied ex- 
tensively [lj] . The TMRG [H method is a powerful 
tool for studying the thermodynamics of one-dimensional 



quantum systems in the thermodynamic limit [14j |. In 
our calculations, the width of the imaginary time slice 
is taken as £=0.1, and the error caused by the Trotter- 
Suzuki decomposition is less than 10 -3 . During the 
TMRG iterations, 60 states are retained, and the trun- 
cation error is less than 10~ 6 . 

Owing to the alternation of the couplings, the TEs of 
the spins coupled by J a and J{ , which are denoted as C a 
and Cf, respectively, are expected to be distinct. Fig- 

(20) 

ure [U[a) shows the temperature dependence of the TE 
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C a at different fields. It is shown that in the absence of 
the field, the intrinsic TE survives, and C a vanishes at 
T"~0.85J a due to thermal fluctuations. In the presence 
of the field, the CT keeps invariant. Although the alter- 
nation is involved, the CT of the intrinsic TE is still a 
fixed point. 

However, the entanglement induced by the field does 
not comply such a rule. For the F couplings, the TE of 
the spins coupled by <// is absent without a field. When 
the applied field closes the gap and increases up to about 
0.9J a , the TE is induced by the field, as shown in Fig. 
[UJb). With further increasing the field, the CT of the 
field-induced TE enhances to reach the maximum at the 
saturation field. A further increase of the field makes the 
field-induced TE vanish as the spins are fully polarized at 
zero temperature. It can be seen that, different from the 
intrinsic TE, the CT of the field-induced TE is dependent 
on the magnetic field. 

Furthermore, the intrinsic TE of the 5=1/2 AF-AF- 
AF-F tetrameric Heisenberg chain [15| is studied using 
the TMRG, which are not presented here. It is shown 
that the CT retains a fixed point, which is also ob- 
served in the trimcrizcd F-F-AF chain 16]. The ob- 
servations suggest that the CT of the intrinsic TE in 
one-dimensional Heisenberg antiferromagnets might be a 
fixed point. 

Next, we treat the AF-F chain within the mean- field 
framework, which may extend the discussions to the 
general 5=1/2 alternating Heisenberg antiferromagnetic 
chains with nearest-neighbor interactions. Following the 
steps in Ref. Ii~7l . we make the Hartree-Fock approxima- 
tion to the Hamiltonian (f2"2"|) after the JW and Fourier 
transforms, and obtain the mean- field Hamiltonian after 
omitting a constant: 

Hhf = J2^ Ja + J t)( d b-\)-h]a{a k 

k 

+ [{J a + Jt){d a - l -)~h]b\bk\ 

+ 5Zt Ja ^ ~ Paby k,2 a{b k + h.c. 

k 

+ Jf(~ ~ Pba)e lk/2 b{a k + h.c], (23) 

where d a ={a} j a. j ) , <4=(&j&i), p &=(&] a j)) Pem=( ffl j+A')> 
which are obtained by self-consistent calculations. Then 
the Bogoliubov transformation is taken to diagonalize 
the above Hamiltonian. Thus, the TE can be calculated 
from the quasiparticle representation. Figure [3] shows the 
mean- field results of the TE C a and Cf. It is shown that 
although the values of the critical fields and CT are not 
accurate, the mean-field results still preserve the features 
of the CT. The intrinsic TE C a vanishes at a common CT, 
while the field-induced TE Cf is dependent on the field. 
As shown in Fig. [3Jb), the CT of Cf enhances with in- 
creasing the field until to the maximum at the saturation 
field, which is analogous to the TMRG result. 
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FIG. 3: (Color online) Temperature dependence of the ther- 
mal entanglement of (a) C a , and (b) Cf for the 5=1/2 AF-F 
chain at various fields obtained by the mean-field theory. 

In this fermion mapping, Z a b=P ab , and X^ b =d a d b — 
Z 2 b , where Z ab and X~T, are the values defined in Eq. ([3]) 
of the spins coupled by J a . Thus, the concurrence C a 
can be expressed by Eq. (|6|) using these quantities. At 
the CT (T c a ), we have 

[M 2 - d a {d b - 1)][M 2 - d a {d b - 1)] = 2\ Pab \ 2 . (24) 

The calculations show that p ab is real, and d a =d b . Thus, 
the above equation can be simplified as 

d a - d 2 a = -V2 Pab - p 2 ab , (25) 

which has the same form as that of the XY chain [Eq. 
(|14|) ] , yielding the following inequality 

+ < \ + <S 2 Vr) 2 (26) 

for the entangled C a . For the field- induced TE Cf, we 
have d a — d a =—^/2p ba —p ba & t the CT, which is satisfied 
at different CTs for different fields. The entangled Ct 

is described by ((5+ 5 2j+1 ) + ^) 2 <i+(5 2 z J ) 2 . Note that 
the thermal quantity witness cannot be written in a form 
as simple as Eq. (|20|) within the present self-consistent 
calculations. 

In summary, we have studied the field dependence of 
the CT of TE in the 5=1/2 spin chains within the ther- 
modynamic limit N-^oo. The concurrence of the TE in 
the spin- 1/2 XY chain is exactly resolved. It is found that 
the CT of the TE is a fixed point. An equation is given to 
determine the CT, which is found to be T c ~0.4843 J and 
smaller than that of the two-qubit system. The thermal 
witness for the entangled state is also proposed. Fur- 
thermore, the TE of an 5=1/2 AF-F chain is studied 
by means of the TMRG method and mean-field treat- 
ment, which indicates that the CT of the intrinsic TE 
of the spins coupled by AF couplings is a fixed point, 
while that of the field-induced TE of the spins coupled 
by F couplings changes with the field. The exact solu- 
tion of the XY chain as well as the mean-field result of 
the AF-F chain indicate that the disappearance of the 
TE is determined by the competition between the spin 
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fluctuations and local magnetic moment at finite tem- 
peratures. The observations suggest that it may be a 
general phenomenon in one-dimensional Hcisenberg an- 
tifcrromagnets that the CT of the intrinsic TE is a fixed 
point independent of the magnetic field. 
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